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Abstract 

Breaking the chiral symmetry, rotation induces a secondary Hopf bifurcation in 
weakly nonlinear hexagon patterns which gives rise to oscillating hexagons. We 
study the stability of the oscillating hexagons using three coupled Ginzburg-Landau 
equations. Close to the bifurcation point we derive reduced equations for the am- 
plitude of the oscillation, coupled to the phase of the underlying hexagons. Within 
these equation we identify two types of long-wave instabilities and study the en- 
suing dynamics using numerical simulations of the three coupled Ginzburg-Landau 
equations. 
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1 Introduction 



Convection has played a key role in the elucidation of the spatio-temporal 
dynamics arising in non-equilibrium pattern forming systems. The interplay 
of well-controlled experiments with analytical and numerical theoretical work 
has contributed to a better understanding of various mechanisms that can lead 
to complex behavior. From a theoretical point of view the effect of rotation on 
roll convection has been particularly interesting because it can lead to spatio- 
temporal chaos immediately above threshold where the small amplitude of 
the pattern allows a simplified treatment. Early work of Kiippers and Lortz 
[1,2] showed that for sufficiently large rotation rate the roll pattern becomes 
unstable to another set of rolls rotated with respect to the initial one. Due 
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to isotropy the new set of rolls is also unstable and persistent dynamics are 
expected. Later Busse and Heikes [3] confirmed experimentally the existence of 
this instability and the persistent dynamics arising from it. They proposed an 
idealized model of three coupled amplitude equations in which the instability 
leads to a heteroclinic cycle connecting three sets of rolls rotated by 120° 
with respect to each other. Recently the Kiippers-Lortz instability and the 
ensuing dynamics have been subject to intensive research, both experimentally 
[4-8] and theoretically [9-15]. It is found that in sufficiently large systems 
the switching between rolls of different orientation looses coherence and the 
pattern breaks up into patches in which the rolls change orientation at different 
times. The shape and size of the patches changes persistently due to the motion 
of the fronts separating them. 

In this paper we are interested in the effect of rotation on hexagonal rather 
than roll (stripe) patterns as they arise in systems with broken up-down sym- 
metry (e.g. non-Boussinesq or surface-tension driven convection with rota- 
tion). The dynamics of strictly periodic hexagon patterns with broken chiral 
symmetry have been investigated in detail by Swift [16] and Soward [17]. They 
found that the heteroclinic orbit of the Busse-Heikes model is replaced by a pe- 
riodic orbit arising from a secondary Hopf bifurcation off the hexagons. Their 
results have been confirmed in numerical simulations of a Swift-Hohenberg- 
type model [18] in which an alternation among the three modes that compose 
the hexagonal pattern is observed. In the following we will call this state 
'oscillating hexagons'. The oscillations can be homogeneous in space or can 
take on the form of traveling waves. Starting from coupled Ginzburg-Landau 
equations, we have previously derived evolution equations for the oscillating 
hexagons that are valid close to the Hopf bifurcation. Within this framework 
the oscillating hexagons were found to support a state of spatio-temporal 
chaos that is characterized by defects, in which the oscillation amplitude van- 
ishes [19]. At the band center the oscillating hexagons and their chaotic state 
are described by the single complex Ginzburg-Landau equation (CGLE). In 
general, however, the oscillation amplitude is coupled to the phases of the un- 
derlying pattern as it is to be expected for a secondary bifurcation. Here we 
study how this coupling affects the stability of the oscillations. In particular, 
it will modify the stability properties of the waves emitted by the spirals in 
the defect chaotic state. We find that the additional coupling leads to new 
long-wave instabilities if the rotation is strong enough or if the wavenumber 
of the hexagons is far away from the band-center.0 

The paper is organized as follows. In the second section we introduce the 



2 For stripe patterns (e.g. convection rolls) the generic equation for a secondary 
Hopf bifurcation is well known [20], and has been studied both theoretically and 
experimentally [21-23]. There, the coupling to the phase of the pattern can delay 
the occurrence of long-wave instabilities of the oscillatory mode. 
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appropriate coupled Ginzburg-Landau equations that describe the hexagon 
patterns and derive the amplitude-phase equations close to the secondary Hopf 
bifurcation. The stability analysis of these equations is addressed in section 
III. In section IV we study numerically the behavior resulting from these 
instabilities. Conclusions are given in section V. Details of the calculations are 
given in two appendices. 



2 Amplitude-Phase equations for Oscillating Hexagons 

We consider small-amplitude hexagon patterns in systems with broken chiral 
symmetry. In order to analyze the possibility of modulational instabilities 
we include spatial derivatives. Due to the strong coupling between modes of 
different orientation, we take the gradients in both directions to be of the 
same order and retain only linear gradient terms (a study of the influence of 
nonlinear gradient terms in the stability of oscillating hexagons is addressed 
in appendix B). After rescaling the amplitude, time, and space we arrive at 
the equations, 

d t A x = ixA x + (m • V) 2 A 1 + A 2 A 3 - A^l 2 (1) 
-[v + i)A x \A 2 \ 2 -{y- 7)^i|A 3 | 2 , 

where the equations for the other two amplitudes are obtained by cyclic per- 
mutation of the indices and \x is a parameter related to the distance from 
threshold. The overbar represents complex conjugation. These equations can 
be obtained from the corresponding physical equations (e.g. Navier- Stokes) 
using a perturbative technique. The broken chiral symmetry manifests itself 
by the cross-coupling coefficients not being equal. Hence 7 is a measure of 
rotation. 

For completeness it should be noted that rotation leads in convection not only 
to a chiral symmetry breaking but also to a (weak) breaking of the translation 
symmetry due to the centrifugal force. In the following we will consider it to 
be negligible. We focus on not too small Prandtl numbers, in which case the 
primary bifurcation is always steady [24,13]. 

Equation (1) admits hexagon solutions Aj = ]Je ,?n J' x+1 ^ with a slightly off- 
critical wavenumber (q^ = + q^, \qj\ <C |q c | = q c ), with 

1 ± ,/l + 4(/i-g 2 )(l + 2z/) 
R= 2(1 + 2.) *-*i + fc + * = °- W 

The stability of this solution to perturbations with the same wave vectors 
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has been studied by several authors [16,17]. Typical results are sketched in 
the bifurcation diagrams shown in Fig. 1. The hexagons appear through a 
saddle-node bifurcation at fj, — /i sn , 



and become unstable via a Hopf bifurcation at /i — with a frequency ujh, 

(2 + 1/) 



+ = 2yfry/{v - if. (4) 



The Hopf bifurcation is supercritical and for \i > stable oscillations in the 
three amplitudes of the hexagonal pattern arise with a phase shift of 2n/3 
between them [16-18], resulting in what we call oscillating hexagons. As \i is 
increased further, eventually a point /i = fihet is reached at which the branch 
of oscillating hexagons ends on the branch corresponding to a mixed-mode 
solution in a global bifurcation involving a heteroclinic connection (see Fig. 
la). Above this point the only stable solution is the roll solution whose stability 
region is bounded below by 

»R = 7 : 7T7 77 + <1 2 - (5) 

(i/ + 7-1)(i/-7-1) 



From Eqs. (4), (5) it is easy to see that, when 7 2 > {y — \){y + + 2), 

the transition to oscillating hexagons occurs at a value of \x for which the rolls 
are still unstable. There is then a parameter regime in which the oscillating 
hexagons are the only stable solution. Furthermore, when [7) > v — 1 = ^ KL 
rolls are never stable and the limit cycle persists for arbitrarily large values of 
fi (see Fig. lb). In the absence of the quadratic terms in Eq. (1) this condition 
corresponds to the Kuppers-Lortz instability of rolls. When the quadratic term 
in Eq. (1) is small (i.e. small non-Boussinesq effects) it can be considered as a 
perturbation of the usual three mode model for rotating roll convection. Far 
above the Hopf bifurcation (fi ^> hh) the periodic orbit is expected to become 
asymmetrical and the resulting state similar to that encountered in the usual 
rotating Rayleigh-Benard convection. 

The stability of steady hexagons with respect to side-band perturbations in 
the presence of rotation has been studied previously [25] . Due to the rotation 
steady and oscillatory, long-wave and short-wave instabilities have been found. 
It turns out that in the presence of rotation the steady hexagons can be stable 
up to the Hopf bifurcation over quite a range of wavenumbers. Thus, the 
oscillating hexagons may be stable near the Hopf bifurcation. 

The focus of the present paper are side-band instabilities of the oscillating 
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Fig. 1. Sketch of a bifurcation diagram for a fixed value of the wavenumber q. a) 
For 7 < 7xl limit cycle of oscillating hexagons ends via heteroclinic cycle and rolls 
become stable, b) For 7 > ^kl oscillating hexagons persist and rolls do not become 
stable. 

hexagons. To address this question analytically we focus on the vicinity of 
the Hopf bifurcation where the oscillation amplitude is small and a weakly 
nonlinear analysis is possible. Since we are dealing with a secondary bifurcation 
the two phases of the underlying hexagons have to be taken into account as 
well. In fact, from a linear stability analysis of Eq. (1) it is easy to see that there 
are four marginal modes at /i — fin- Two correspond to the Hopf bifurcation 
and can be described by a complex amplitude Ti. The other two correspond 
to a two-dimensional phase vector 0, related to translations in the x- and 
y-directions [26,27]. It can be written as a combination of the phases of the 
three modes: = ((f) x ,(j)y) = (—02 — 03) (02 ~ 03)/v / 3), satisfying the locking 
condition 0i + 2 + 03 = $. The global phase $ is a fast variable that relaxes 
rapidly to its stationary values $ = 0, ir (up or down hexagons). To study the 
nonlinear behavior of the oscillating hexagons, the amplitudes Ai are expanded 
as: 

A n = (R H + e 2 ™ /3 ^e !Wff * + c.c. + 0(e))e l ^ + ^ , (6) 



where e is a small parameter related to the distance from the bifurcation line. 
The amplitude of the steady hexagons at the bifurcation point, Rh = 
is independent of fi and q. Eliminating the fast variables, at order e 3 / 2 (see 
appendix A) we arrive at an equation for the amplitude of the oscillation 7i, 
coupled to the phase vector of the underlying hexagonal pattern, 

d T n = e5{H + iV 2 n- 5 2 1-LV -0- pH\H\ 2 , (7) 
( 9 T = J D ± V 2 + ^!|V(V-0) + J D Xl (e 2 x V 2 0) (8) 
+D X2 (e z xV)(y ■ 0) + aV\H\ 2 + A(e, x V)\H\ 2 
-i[3 2 (HVH - HVH) + ir)[H(e z x V)H - H(e z x V)W], 
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with e = jj, — hh and 



v = 3R H (l + 2R H ), 
2R H Huju 



Si 

P = 
D ± 

a = 



v v 
1 ?~ 2 



S 2 = q5 u 
3q 2 R H iq 2 9R\ 



2u? H 



2 9R 2 H + u 2 H u H 9R 2 H 
8(3R H + 1) Mu H {\ + 4R H ) 32i 



R H v 



3u H 



1 1 

= p ^11 = 2 

9i% + c4 " 
R H (9R 2 H + u 2 H ) 



_2y 

2g(l + 6i2 H ) 



R H v 

, #2 = A, ^ 



X2 



18g 



0, 



9R 2 H + u 2 H 



(9) 
(10) 

(11) 

(12) 

(13) 

(14) 

(15) 



It is worth pointing out that the phase-amplitude equations (7,8) can be de- 
duced by means of symmetry arguments alone and are, therefore, generic to 
this order in e. In fact, they could be derived directly from the fluid equations 
without the use of the Ginzburg-Landau equations (1). Thus, keeping higher 
order terms in (1) would change the values of the coefficients, but not their 
form. 

When deriving (7,8) using symmetry arguments, one interesting aspect has 
to be taken into account. In most secondary bifurcations, the oscillating am- 
plitude is either even or odd under reflection symmetry (see, for instance, 
[20]). In our case, however, the field H transforms under reflection x — > — x 
as Ji — * 7i. The temporal phase of this complex amplitude is therefore a 
pseudo-scalar, changing sign under reflection. This is because this phase is 
related to the oscillating frequency which, in turn, depends linearly on the 
chiral symmetry breaking coefficient {ujh = 2\/3R 2 H 'y). This implies that, in 
Eq. (8), the term (HVH — HVH) breaks the chiral symmetry, but not the 
term [H(e z x V)7i — 7i(e z x V)7i], as one could have naively expected. Looking 
at the values for the coefficients we see that, in fact, /3 2 changes sign while rj 
is invariant under uj h — > — uj h . 

It is interesting to note that at the band-center (q = 0) the system (7,8) 
decouples. In this case the usual CGLE for the amplitude of the oscillation is 
recovered, which after rescaling the amplitude, time, and space can be written 
as [28]: 

d t H = H + (l + ib x )V 2 H - (6 3 - i sign{u H ))H\H\ 2 , (16) 
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with 



6 2(i% + 2^)g 2 

01 £ r (2qm H -Rl-ujj 1 )uj H ' {U> 

, Pr _ 2\u H \R H (3R H + 1) ( 

3 " + 4i2 H )+8i&(l + 2i2 H )' 1 J 

where the sub-indices r and i indicate real and imaginary part, respectively. 
Furthermore, at the band-center b\ — 0. 

The CGLE has been studied extensively. It possesses an extraordinary variety 
of solutions, including a phase chaotic state, defect chaos, a frozen vortex state 
and stable plane waves [28-35]. For the case considered here the values of the 
parameters b\ and 63 are always such that the system is in a regime in which 
stable plane waves coexist with defect chaos. In the present context it has 
to be emphasized that the complex amplitude of the CGLE represents the 
amplitude of oscillation of the hexagon pattern. Thus, a solution of Eq. (7) 
with spatially uniform amplitude Ti = He tnt , 



corresponds to a state in which all three amplitudes oscillate in time phase- 
shifted with respect to each other but the phase of each amplitude is constant 
in space, as illustrated in Fig. 2. In a traveling wave solution (TW), 

H = He int+ikx , (20) 



with 



H 



Ie8i r — (l r k 2 

Pr 



n = e(5u-^5 lr )-^ t -^ r )k 2 , k=\k\, (21) 



the phase of each amplitude is space dependent and in different parts of the 
system roll-like hexagons with different orientation are dominant at any given 
time. This is shown in Fig. 3. Note that such a state has two different wavenum- 
bers: that of the underlying regular hexagon pattern and that of the wave 
modulating the oscillation amplitude. At the band-center, i.e. for hexagons 
that have the critical wavenumber, the modulation of the oscillation ampli- 
tude does not affect the phase of the underlying hexagons. 

Away from the band-center (q 7^ 0) the complex amplitude Ti and the phase 
vector become coupled. For the traveling-wave state of the oscillation am- 
plitude (20) this implies that even the underlying hexagon pattern drifts with 
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Fig. 2. Oscillating hexagons (/x = 4.6, 1/ = 2, 7 = 0.5, L = 200, k = 0, q = 0, 
g c = 0.44). There is one period of oscillation between the first and last snapshots. 
The contour lines are taken at ip = J2j=i Aje 1 ^'* = —1, 0, 1. 
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Fig. 3. Traveling hexagons (ji = 4.6, v = 2, 7 = 0.5, L = 200, fc s = 0.0628, 
k y = 0, q = 0.125, q c = 0.31). There is one period of oscillation between 
the fist and last snapshots. The modulation is traveling with phase velocity: 
v p h = (ooh + ~ 27.27. The underlying hexagons are traveling with a speed: 

v = 2H 2 ^J ft 2 + V 2 /Q — 0.00035, too small to be observed in the figures. 
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a velocity that depends on the wave vector of the modulation. Specifically, the 
phase is given by = ut, with 

C3=-2H 2 [p 2 k-r ] (e z xk)}, (22) 



which implies a drift of the hexagons with a speed v = 2H 2 ^J(3 2 + r] 2 /q, at 
an angle 9 = arctan(— 77/ /? 2 ) with respect to the wave vector of the traveling 
wave. Substituting the values of the coefficients (9-15) one obtains for the 
angle 6 = arctan(3i?if/u;#). Thus, for u>h — > the angle becomes 6 = ±7r/2, 
i.e. the drift is perpendicular to the wave vector of the modulation. When 
u H — > 00, on the other hand, they drift in the parallel direction (0 — > 0, n). 
The speed v is given by 




The traveling waves exist above the curve eSi r — £ r k 2 = 0, up to the global 
bifurcation at fj, — \ihet- However they can be unstable to side-band pertur- 
bations. In the following we study their stability properties. In particular, we 
will focus on the effect of the coupling of the oscillation amplitude Ti to the 
phase of the underlying hexagons. 



3 Linear stability analysis 

To consider the linear stability of the oscillating hexagons, we perturb the 
traveling-wave state (20) as 

H = (H + h)e l(nt+k - x+ v\ <P = ut + 0i. (24) 

For simplicity we write 0i in the following as <fi. Substituting (24) in Eqs. (7,8) 
we obtain the linear equations for the perturbations: 

d T h = -2£ r Hk ■ Vy? - 2&k • Vh + £ r V 2 /i - &if VV - 5 2r HV ■ <f (25) 
-2 Pr H 2 h, 

d T cp = -2&k • Vy? + %k • Vh + f r V V + Tf y2/l - 5 2i v • $ ~ 2piHh(26) 

H H 

d T $= D ± V 2 $ + £>||V(V • 0) + £> Xl V 2 (e 2 x 0) + D X2 (e 2 x V)(V-0)(27) 

+2aHVh + 2p 1 H(e z x V)h - 2(3 2 H 2 V^ + 2r]H 2 (e z x V)<p 

-4(3 2 Hhk + 4r]Hh(e z x k). 
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This leads to a 4 x 4 linear eigenvalue problem, which must be solved numer- 
ically. In the long-wave limit the perturbation in the amplitude h becomes 
slaved to the gradients of the phases ip and <f>. The resulting 3x3 system is, 
however, still rather involved. A substantial simplification occurs in the case 
k = 0, i.e. for homogeneously oscillating hexagons. In order to gain physical 
insight we will consider this case first. 



3.1 Homogeneous Oscillations 



For homogeneously oscillating hexagons u = and the perturbation equations 
(25,26,27) reduce to 



d T h = i r V 2 h - £HV V - S 2r HV ■ 4> - 2p r H 2 h, (28) 

d T V = + %V 2 h - 5 2i V ■ - 2 Pl Hh, (29) 

d T $= D ± V 2 $ + £>||V(V • $) + £ Xl V 2 (e 2 x $) + D x , 2 {e z x V)(V-0) 

+2aHVh + 2(3 l H(e z x V)h - 2(3 2 H 2 V^ + 2r]H 2 (e z x V)</?. (30) 



In the limit of long-wave perturbations the amplitude h can be eliminated 
adiabatically 

h ^ ^Jpl~^ H ^ - 5 *rHV ■ # (31) 



In this manner we arrive at a system for the three phases, two corresponding 
to the spatial translations in the plane, the other to a temporal shift, 



d T $ = D L V 2 $+ (d„ - ^y- j V(V • 0) + D Xl V 2 (e 2 x 
+ (d X2 - S^l) (e z x V)(V • 0) - ^V 2 (Vy?) 

V Pr J Pr 



-V 2 (e 2 x V)v? - 2(3 2 H 2 Vif + 2??iJ 2 (e 2 x V)y?, (32) 



Pr 

dr¥> = Ur + V 2 y? + UA ~ S 2t ) V • 0. (33) 



At the band-center, q = 0, Eqs. (7,8) decouple since /5i = (5 2 = S 2 = a = r\ = 
(cf. (9-15)). It is easy to show that in that case the eigenvalues are 
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^1,2 = "2 



2D ± + D l{ ± ^Dg - 4D Xl (D Xl + L> X2 )j Q 2 , (34) 

= - Ur + Q 2 - (35) 



We obtain therefore the usual expression for the phase instabilities of the 
underlying hexagons [25] and the Benjamin- Feir instability of the oscillations 
[36]. The actual values of these eigenvalues, when q — 0, are o\ = —Q 2 /4, 
a 2 = — 3Q 2 /4 and <7 3 = —Q 2 /2. The system is therefore always stable at the 
band-center. 

Away from the band-center (q ^ 0) the system is no longer decoupled. At 
leading order in the long- wave expansion (32), (33) we have then 



d T (p = -2(3 2 H 2 V V + 2 V H 2 {e z x V)y?, (36) 
& r <p=(62rj-5 2t y jV-$. (37) 

These two equations can be combined into a single second-order equation for 

I 

dU= "2 UA - 5 2 ^j H 2 (fcV - V (e z x V))(V • 0). (38) 



Writing in normal modes, <f> x = 0°e iQ ' x+at , <p y = 0° e iQ - x+<T *, <p = (p°e iQ - M , 
Q = |Q|, we arrive at the dispersion relation 



<7i = 0, (39) 



a 2 , 3 = ±p (5 2rPi /p r - hi) fo HQ, (40) 
indicating the possibility of two different instabilities. 

The eigenvalue o\ corresponds to the divergence-free part of <fi and does not 
involve ip (cf. Eq. (37)). It is marginal at this order. The eigenmodes associated 
with <r 2 ,3 involve both phases and ip. For 7 2 > 2/(?>R 2 H ), a 2 is always positive 
(except at the band-center, where (3 2 = 0). Hence there exists a critical value 
for the rotation \~f c \ = ^2/3/R H above which the system is unstable for Q — > 
0. When |7| < |7 C | the eigenvalues a 2 $ are purely imaginary and the stability 
is determined at next order. 

At quadratic order in the perturbation wavenumber Q one obtains 
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Substituting the values of the coefficients from (9-15) in <j\ we see that it 
becomes positive provided q 2 > 7 2 -R#. In Fig. 4a the long-wave stability limits 
of the oscillating hexagons are shown for v = 2 and 7 = 0.5. These results 
agree with those obtained solving the full 4x4 dispersion relation. Hence, 
for any given value of the rotation rate there is a value of the wavenumber 
above which the system becomes unstable. The range of stable wavenumbers 
decreases as the rotation rate is decreased. In fact, it vanishes as 7 — ► 0. In this 
limit also the range in e over which the oscillating hexagons exist vanishes. 
Also shown in Fig. 4a are the instabilities of steady hexagons, below the Hopf 
bifurcation. The solid line represents the long-wave results, while the circles 
are the results obtained solving the 6x6 dispersion relation associated with 
Eq. (1) [25]. The dash-dotted line represents the line above which the rolls 
become stable (cf. Eq. (5)). 

As 7 is increased the range of wavenumbers that are stable with respect to the 
diffusive mode {a\ < 0) increases. However, for 7 > j c , cr 2 becomes positive 
for Q — > 0. In finite systems the term quadratic in Q may not be negligible. In 
fact, inserting the values from (9-15) into Eq. (42) shows that the g 2 -term in 
(72,3 is always stabilizing when q <C 1 (cr 2 = 03 = — 5Q 2 /8, when q = 0), and 
for q ~ 0(1) it is only destabilizing when 7, v S> 1. For q = 1 typical values 
of 7, v for which this happens are v ~ 30, 7 ~ 70. Therefore, in finite systems, 
in which Q m i n = 2tt/L cannot be arbitrarily small, there is always a region 
close to the band-center that is stable, even when I7I > |7 C |. The stability 
limit cr 2 = is given by an expression of the form if 2 ~ fi — fin = f(q) / (Lq 2 ), 
where f(q) appears due to the g-dependence of the second term in (42). This 
situation is shown in Fig. 4b, where the different symbols correspond to several 
values of the system size. For smaller systems the stable region increases. The 
shape of the stability limits suggest that f(q) does not depend strongly on q. 
In Fig. 4b the instability corresponding to o\ > has moved to higher values 
oiq (q = ±0.9).Q 



3 The form of the eigenvalues (12,3 is similar to that encountered in the case of a 
secondary Hopf bifurcation off a roll pattern [21,22]. However, in contrast to the 
one-dimensional case there is a third eigenvalue, a\, because the phase has two 
components. Similar expressions for the eigenvalues are to be expected for other 
two-dimensional patterns, e.g. squares, undergoing a Hopf bifurcation. 
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Fig. 4. Stability limits of oscillating hexagons [y = 2, j c = 0.816). Lines indicate 
saddle-node for hexagons (dotted), Hopf bifurcation (dashed) and stability limits 
for rolls against hexagons (dot-dashed), a) 7 = 0.5 < j c . b) 7 = 0.9 > j c . Symbols 
indicate different system sizes, L=250 (circles), 500 (squares) and 1000 (diamonds). 
The stable region close to the band-center shrinks as the system size is increased. 

It should be noted that taking the limit q — * in (41), (42) does not give the 
same results as in (34), (35). In fact, in the limit q — > 0, the three eigenvalues 
in (41), (42) become a 1 = -Q 2 /4, a 2 = a 3 = -5Q 2 /8, while Eqs. (34), (35) 
yield o 2 = —3Q 2 /A and 03 = —Q 2 /2. This difference is due to the fact that, 
in order to obtain (41) and (42) we assume q to be 0(1), and expand in terms 
of Q, while (34) and (35) are obtained by taking the limit q — > first. If we 
consider both to be small and of the same order, q ~ Q 1, then: 



<Ji = ~D±Q 2 
^ = -\{D ± + D\\+i r )Q 2 



Pi 



(43) 
(44) 



1/2 



Q. 



When q = we obtain o\ = -Q 2 /4:, a 2 = -3Q 2 /A and 0-3 = -Q 2 /2, while 
considering Q q expressions (39), (40) are recovered at leading order. 



3.2 Stability of Traveling Waves 



For the traveling waves the expressions become quite complicated. We there- 
fore present the analytical results from the long-wave analysis only up to linear 
order in the gradients. As in the case of homogeneous oscillations, in the long- 
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wave limit h becomes slaved (cf. (25)), 

1 



h ~ 



2p r H 2 



[-2{; r H{k-V)y-5 2r H{V •</>)]. 



(45) 



At the band-center and decouple and we recover expression (34) for the 
eigenvalues 0\^ 2 associated with <fi. The eigenvalue a 3 describes the long-wave 
behavior of the single CGLE. It is given by 



a 3 =iv g -Q + 2(1 + ~ & + ^)Q 2 , 

Pr to lr ~ S>K< Pr 



(46) 



and yields the usual Eckhaus stability limit for waves. Here we have introduced 
the group velocity 



_dtt pi 



(47) 



The expression e = p — pu = £, r k 2 /8i r gives the neutral surface for the ap- 
pearance of traveling waves of wavenumber k. 

Away from the band-center Eqs. (25,26,27) reduce at leading order in the 
long-wave expansion to 



<9 T = ^[/3 2 k - V (e z x k)](k • V)y? + —Ifck - V (e z x k)](V • $) (48) 

Pr Pr 

-2H 2 [p 2 V<p-rj(e z xV)<p], 



d TV = -2(& - %)k .V<p+[5 2 A-5 2i )V-l 

Pr V Pr J 



(49) 



This yields the eigenvalues 



(71 =0 

*2,s = £(v, + ^v fc ).Q± 4 



v s + — v h I • Q 

Pr J 



1/2 



+-(5 2r v g • Q - Cr-rk • Q)(v ff • Q) + 2H 2 f5 2 TQ 2 

Pr 



(50) 
(51) 



where we have defined T = 5 2r pi/ p r — 5 2 « and v ft = —uo/(2H 2 ) (cf. (22)). 

Since the analytical expressions for the eigenvalues are quite complicated we 
solve Eqs. (25,26,27) directly without the additional long-wave approximation 
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Fig. 5. Stability limits of traveling waves [y = 2, j c = 0.816). a) 7 = 0.3 < 7 C . 
Shown are the neutral surface and the stability limit. Traveling waves are stable 
above the second surface, b) 7 = 0.9 > j c . Dashed-dotted line: stability limit of 
rolls against hexagons. 

(45) for various values of the wavenumber k = (k, 0) of the traveling waves. 
Fig. 5a shows the neutral surface of traveling waves and their stability limit 
as a function of the hexagon wavenumber q. Due to the reflection symmetries 
q — > — q and k — > —k only one quadrant is shown. For clarity the stability 
surface has been capped at fi = 6. For k = the stability limit does not 
depend on jj, and the stability surface is vertical (cf. Fig. 4a). As k is increased 
the stability surface becomes smoother. For k 7^ 0, the waves are unstable 
at onset and become stable above the second surface. Fig. 5b shows cross- 
sections of the stability surface for 7 > 7 C , for a system of size L = 250. When 
k = we recover the results from Fig. 4b. For k 7^ 0, but small (cf. k — 0.1 
in Fig. 5b) the traveling waves are unstable at onset but they become stable 
for larger values of /i. When \q\ is large they can become unstable again as 
fi is increased. For larger values of k, the latter two stability lines merge and 
the stability region becomes bounded in q. At this point the stability limits in 
Figs. 5a and 5b look similar. 



4 Numerical simulations 



In order to study the nonlinear behavior arising from the instabilities, we 
have performed numerical simulations of Eqs. (1) and (7), (8). A Runge-Kutta 
method with an integrating factor that computes the linear derivative terms 
exactly has been used. Derivatives were computed in Fourier space, using 
a two-dimensional fast Fourier transform (FFT). The numerical simulations 
were done in a rectangular box of aspect ratio 2/y/3 with periodic boundary 
conditions. This aspect ratio was used to allow for regular hexagonal patterns. 
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Fig. 6. State resulting from the long-wave instability a\ > (y = 2, 7 = 0.3, g = 0.5, 
/x — /i/f = 0.1, L = 250). a) The dashed and solid lines mark the zero contour lines 
of the real and imaginary part of the oscillation amplitude Ti. b) Reconstruction of 
the hexagonal pattern from a), with q c = 0.4 

We investigate the stability of oscillating hexagons simulating both the orig- 
inal amplitude equations (1) and the reduced amplitude-phase equations (7), 
(8). To that end, we start with homogeneously oscillating hexagons as given 
by Eq. (19) (or Eq. (6)) and add weak noise. For values of /x close to fi H , the 
growth rates obtained from (1) and (7), (8) agree with each other and with 
the results from the linear stability analysis. Both the long-wave instabilities 
coming from o\ > (41) and 0-2,3 > (42) lead to qualitatively similar be- 
havior. The perturbations grow until they reach a saturation suggesting that 
the bifurcations are supercritical. Since the perturbations involve spatial mod- 
ulations of the oscillation amplitude the hexagons begin to travel. However, 
the modulation wavevector k varies in space and induces drift velocities that 
are different in magnitude and direction at different locations in the system, 
implying a shear of the pattern. This results in deformed hexagon patterns as 
shown in Fig. 6. 



In addition to the long-wave instabilities a short-wave instability appears for 
larger values of the hexagon wavenumber q. This instability is induced by 
the short-wave instability of the steady hexagons [25] and cannot be studied 
with the amplitude-phase equations (7), (8). As /i is increased above fin, the 
oscillating hexagons arise through a Hopf bifurcation off the steady hexagons. 
Since that bifurcation occurs at zero wavenumber it affects the long-wave 
properties of the system and we expect that the long-wave stability limits of 
the steady and of the oscillating hexagons differ qualitatively. The effect of the 
Hopf bifurcation on short-wave instabilities, on the other hand, should vanish 
as the amplitude of the oscillations goes to 0. Thus, we expect a continuous 
transition from the short-wave instability of the steady to that of the oscillating 
hexagons as the line fi = fin is crossed. In order to check this, we numerically 
determine the short-wave instability of the oscillating hexagons and compare it 
with the stability results for the steady hexagons, as obtained from solving the 
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Fig. 7. Short-wave stability of oscillating hexagons (stars) and steady hexagons 
(circles) for the same parameters as in Fig. 4a [y = 2, 7 = 0.5). Solid line: Hopf 
bifurcation of steady hexagons. 

sixth-order dispersion relation associated with Eq. (1) [25]. The transition is 
indeed continuous, as expected. The stability region of the oscillating hexagons 
turns out to be reduced as compared to that of the steady hexagons when \i 
is increased (see Fig. 7). 

If the noise added to the oscillating hexagons is sufficiently large, different dy- 
namics may arise even in the parameter range in which the oscillating hexagons 
are linearly stable. As indicated earlier, within the framework of the three 
coupled Ginzburg-Landau equations (1) one obtains for the parameters in the 
CGLE (16) values for which a persistent chaotic state exists while the plane 
waves are linearly stable. A detailed study of the chaotic state, which is char- 
acterized by the creation and annihilation of defects, comparing its description 
using the coupled Ginzburg-Landau equations (1), the amplitude-phase equa- 
tions (7,8), and the single CGLE (16) has been presented elsewhere [19]. One 
of the main results of that study is the observation that this system is one of 
the few in which the defect chaos of the CGLE should be accessible experi- 
mentally. 



5 Conclusions 

In this paper we have studied the effect of a breaking of the chiral sym- 
metry on systems that exhibit hexagon patterns. Classic examples of such 
systems are non-Boussinesq and surface-tension-driven convection with rota- 
tion. We have focused on the dynamics of the oscillating hexagons that arise 
in a Hopf bifurcation due to the rotation. In the vicinity of this secondary 
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Hopf bifurcation the oscillating hexagons are described by a single complex 
Ginzburg-Landau equation (CGLE) coupled to the two phases of the underly- 
ing hexagons. The resulting amplitude-phase equations have certain similari- 
ties with those describing the secondary Hopf bifurcation observed in rectangle 
patterns in electro-convection in nematics [23]. 

Like the CGLE, the amplitude-phase equations support homogeneously oscil- 
lating solutions as well as traveling waves. In the latter, the coupling to the 
phases induces a drift of the hexagons in a direction that is typically oblique to 
the propagation direction of the traveling waves. The stability analysis of the 
oscillating and the traveling hexagons reveal two types of long-wave instabili- 
ties, one occurring when the hexagons are far away from the band-center, the 
other for high enough rotation rate. Even in this latter case, in finite systems 
there is always a stable region close to the band-center, its size depending on 
the size of the system. In both cases, the instabilities appear to be supercriti- 
cal, giving rise to a spatially modulated oscillating hexagonal pattern. 

Although there is always a region in which the homogeneously oscillating 
hexagons are linearly stable within the three coupled Ginzburg-Landau equa- 
tions, they are in fact only meta-stable. Sufficiently large perturbations induce 
a transition to a state of defect chaos described by the CGLE [19]. There is 
always bistability between the chaotic state and the ordered oscillations. In 
that respect the chaotic state resembles spiral-defect chaos as it is observed in 
Rayleigh-Benard convection at low Prandtl numbers [37]. As in that system 
the ordered states can presumably only be obtained by carefully controlled 
initial conditions. Rotating hexagon convection appears to be the first system 
in which the defect chaotic regime of the CGLE should be accessible experi- 
mentally. 

From numerical simulations of the amplitude equations (1) it has been shown 
that a transition between defect-chaos and a frozen vortex state occurs for 
wavenumbers far away from the band-center [19]. This transition happens 
when the asymptotic plane waves emitted by the spirals become absolutely 
unstable. Therefore, the stability results for the traveling wave solution are 
relevant in order to determine the range of existence of the chaotic state found 
in [19]. A complete treatment of this point would involve the determination of 
the asymptotic wavenumber koo for the waves emitted by the spiral solutions 
of the coupled amplitude-phase equations. 

In the present work, we have studied the system close to the bifurcation point, 
where analytical results can be obtained. However, far away from the Hopf bi- 
furcation, a number of interesting effects are expected. One is related to the 
strength of the chiral symmetry breaking. To leading order in the perturba- 
tion expansion, the usual single CGLE is obtained, which is chirally symmetric. 
The breaking of the chiral symmetry appears only through the coupling with 
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the phase, and vanishes at the band-center. For larger values of the oscillation 
amplitude, higher-order terms breaking the chiral symmetry of the CGLE are 
expected. This leads to an asymmetry between defects with opposite topo- 
logical charge [38], with one type of spirals becoming dominant in the frozen 
vortex state. The asymptotic wavenumber, as well as the onset of absolute in- 
stability, becomes different for positively and negatively charged spirals, and 
the transition to the defect chaotic regime is different than in the chirally sym- 
metric case [38]. We expect that this effect of the chiral symmetry breaking 
can be observed with Eqs. (1), as the control parameter is increased above the 
Hopf bifurcation. 

Another open question is the relation of the limit cycle corresponding to the 
oscillating hexagons with the heteroclinic orbit arising in the Kiippers-Lortz 
instability of rolls. As the limit cycle approaches the mixed-mode solutions, the 
harmonic oscillations are transformed into a switching between the three roll 
modes making up the hexagons similar to the dynamics arising in the Kiippers- 
Lortz regime. This suggest a connection between the domain chaos found in 
these systems and the regular and disordered states discussed in this paper. 
However, it is worth emphasizing that in the defect chaos regime described 
in [19], the orientation of the hexagons is well defined and what is spatially 
chaotic is the modulation of the amplitudes that compose the hexagon pat- 
tern. In the Kiippers-Lortz regime, on the other hand, patches of rolls with 
arbitrary orientation are possible, resulting in a state with an isotropic Fourier 
spectrum. Thus, a quantitative comparison between both states is not possible 
with Ginzburg-Landau equations such as (1) and generalized Swift-Hohenberg 
models must be considered [39]. 
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A Derivation of the Amplitude-Phase Equations 

At the Hopf bifurcation, the hexagon solution (2) becomes Rh = \j[y — 1). 
We will consider perturbations around this solution, both in amplitude and 
phase: 

A t = (i^ + r,)e^ x+ ^+*, (A.l) 
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where $ = 0i + 2 + 03 is the global phase of the hexagons and the three 
phases 4>i can be written as 

01 =&x, 

02 = -0^/2 + V3(j) y /2, 

03 = -0^/2 - V3<j) y /2, 

where we have defined the phase vector = (<p x ,(p y ), with X and y related 
to translations of the pattern in the x- and y-directions. 

For the perturbations of the modulus there are three eigenvalues, one real, 
U\ = —2R 2 (1 + 2v) + R, corresponding to an eigenvector with r\ = r 2 = r 3 
and a complex conjugate pair, <t 2 ,3 = — 2i? 2 (l — v) — 2R ± 2y / 3i? 2 72, whose 
real part vanishes at /i — fin- The corresponding eigenvector satisfies r 3 = 

e 2 -/ 3 r 2 = e 47ri / 3 n. 

Taking this into account, we consider the expansion: 





= er + 


{-\ + i^) WlHe^ + e(Wi + ^i 2 e 


- 2 ^'))+c.c. 


r 2 


= er + 


(-^ - i^j-)(V~ene^ + e(Wi + W 12 e 


- 2 ^)) + c.c. 


r 3 


= er + 


XyftHe™ + e(Wi + W 12 e- 2 ^)) + c.c." 


J 


<t>x 


= Ve0x 


+ e(0 :rl e^ + C . C .), 




02/ 


= v /e"0j / + e(0 ?/1 e^ + c.c.), 




$ 


= e($ie 


^ + c.c). 





We also assume that the resulting state evolves on long time and space scales, 
specifically: d t ~ 0(e), V ~ (^(e 1 / 2 ). Substituting the former expressions into 
Eq. (1), at (^(e 1 / 2 ) the linear problem is recovered, giving the value for the 
critical frequency: ujh = —2^/3Rj I 'j. 

At 0(e) an algebraic relation between the slaved and the marginal modes is 
obtained: 



0xi = ^|-(v / 3 + i)(d x H - id y H), (A.2) 
^ = -^(^ + - id y H), (A.3) 

$l = 9P ^g g ; ^ - + (A - 4) 
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1 



n 



3R H {1 + 2R H 
qR H 



HR H - (1 + 6R H )\H\ 2 - qR H V 



10 



— (VSi + 1) [d x (j)y + dy(j) X + i(d X (f) X - dy(j)y)} , 

H 



« a = 5(l + -l«=. 

6 V o>#, 



(A.5) 
(A.6) 
(A.7) 



At (9(e 3 / 2 ) we obtain a solvability condition for H and <p: 



d t H = fiH+ -V 2 H-qHV 



1 



i2? 



H 

R H q 



UJl 



3R„(l + 2R H )- — 



\n\ 2 n + (8-^-)nn 12 



R 



ii 



1 + >/3i) [fte^xi - dy(pyi + i(dx<Pyi + d y (j) xl )} 



+ ^Y-(V3i + l)(d X <t> 1 -idy<t> 1 ), 

d t <f> x = ]v 2 <p x + \d x {V • 0) + ^-d x r 
4 Z Kh 



q 



2R H 



4R H 
1 



2R 



H 



(VSi - l)(d x Hio - idyTLw) + c.c. 



- V3)H((f) x i + i(f) y i) + c.c. 



(1 + V3i)(3R H - ui)H$i + c.c. 



dt<Py=\v 2 <t ) y+ 1 -d y (V-$) + ^-dyr 



2R H 

uj 



[( VS + i) (d x H w - idyHva) + c.c. 
[(1 + y/3i)H((j> x i + i<f> y i) + c.c. 
[{V3 - i)(3R H - cji)H$i + c.c. 



(A.8) 



(A.9) 



(A.10) 



Substituting the expressions of the slaved modes into the former equations we 
obtain the amplitude-phase equations (7), (8). 



B Nonlinear Gradient Terms 



If we retain in the Ginzburg-Landau equations (1) the nonlinear gradient terms 
that express the dependence of the quadratic coupling term on the hexagon 
wavenumber [25] we obtain 
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d t A ± = nA x + (m • VfA 1 + A 2 A 3 - A ± \ A ± \ 2 

-(u + ^A^-iu-^AMsl 2 

+i{ai + a)A 2 (n 3 • V)A 3 + i( ai - a)A 3 (n 2 ■ V)A 2 

+ia 2 (A 2 (t 3 ■ V)A 3 - A 3 {t 2 ■ W)A 2 ) . (B.l) 
For the amplitude-phase equations (7), (8) we obtain then the coefficients: 



a = 1 + 2qai, 
B - a 

V — 1 

a 2 (2 + u) 



° (^-l) 2 
v = 3R H (a + 2R H ), 

2aRn 2iuj h 

Oi = 



+ q 2 = Rh(3Rh + a) + q 2 , 



V V 

2aRnq AR 2 h oii{3Rh + oi) 2iu>H(q — Rh&i) 

. - - - -. — — .- " " ■ 



V V V 

1 R H (R H a 1 + q) V3 „ R H r y/3R H qa 



i(R H ai + g) 
9R 2 H a 2 + uj 2 H 
iR 



Rh 

uhQ + — [wh{oli - V3a 2 ) - 3V3R H aa] 



iq 



H 



H 



[oi\ + V3a 2 ) 2 + 3a 2 



Sa(3R H + a) Aiuj H (a + AR H ) 32ia 2 



R H v 



3u 



H 



Dh = 



1 \f3Rnqct 

4 u H 

1 RnOii — q 

2 v 



Rh(&i - V3a 2 ) - 2q 



D X1 = 



Rjj 



(«i + V3a 2 ) 2 + 3a 2 ] , 



D X2 = v ^ Rh ° : (R h01i - q), 



a = -^T [2q + Ru{ ^ a2 ~ ai)] + RU9RW + . 2 ;) ' 
6u H a(q + R H ot\) . \[3a 

A = n /nn2.9 , .2 \ H ^ + 6 ^)' 



(32 = 



R H (9R 2 H a 2 + lu 2 h 
Qoj H Ci{q + RhOli) \f3< 



R H (9R 2 H a 2 + u 2 H ) R H ' 
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_ [9R 2 H a 2 (2q - R H {u x + y/3a 2 )) - R H u} 2 H (^ a i + v 7 ^)] 
V ~ R 2 H (9R 2 H a 2 + u 2 H ) 

B.l Comparison with the CGLE 



After rescaling we obtain the values for the coefficients of the CGLE (16). The 
coefficient 63 now becomes: 



2\ui H \R H a{a + 3Rh) 
uo 2 H (a + 4R H ) + 8R 2 H a 2 {a + 2R H ) ' 



°3 = ..2 ,„ ~ A D N I OD2 n 2(~. 1 o D \ ' ( B - 2 ) 



For q 7^ the maximum value of 63 is: 

bT x = a + 3RH (B.3) 

^8(a + 2R H )(a + 4R H ) 

and the limits for small and large Rh are the same as in (18), even if now 
63 depends on q (through a — 1 + 2qct\). At the band-center 63 becomes the 
same as (18). 

The expression for b\ is quite involved. At the band-center it reduces to: 
_ {9R 2 H + ujj)[( ai + V3a 2 ) 2 + 3a 2 ] + 2uj H a 1 [uj H {V3a2 

1 = / = 7 =— •( / 

2u H {9R 2 H + uj 2 h - ai-R^[3i?//(ai - V3a 2 ) + v3aw ff ]) 

An important change occurs with respect to the decoupling of the amplitude- 
phase equation. For «! ^ the amplitude and the phases do not decouple 
anymore. It would require <5 2r = 0, <5 2i = 0. From the latter we obtain q = 
Rh&i, while the former implies 



-(1 + 4R H al) ± J (I - 4i? H «i) 2 - 48R 2 H a 2 



4ai 

The only real solution for both conditions is a,\ = 0. 
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